Quantumness in decoherent quantum walk using measurement-induced disturbance 
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The classicalization of a decoherent discrete-time quantum walk on a line or an n-cycle can be 
demonstrated in various ways that do not necessarily provide a geometry-independent description. 
For example, the position probability distribution becomes increasingly Gaussian, with a concomi- 
tant fall in the standard deviation, in the former case, but not in the latter. As another example, 
each step of the quantum walk on a line may be subjected to an arbitrary phase gate, without 
affecting the position probability distribution, no matter whether the walk is noiseless or noisy. 
^—^ , This symmetry, which is absent in the case of noiseless cyclic walk, but is restored in the presence 

of sufficient noise, serves as an indicator of classicalization, but only in the cyclic case. Here we 
show that the degree of quantum correlations between the coin and position degrees of freedom, 
quantified by a measure based on the disturbance induced by local measurements [Luo, Phys. Rev. 
A 77, 022301 (2008)], provides a suitable measure of classicalization across both type of walks. 
Applying this measure to compare the two walks, we find that cyclic quantum walks tend to classi- 
cs] | calize faster than quantum walks on a line because of more efficient phase randomization due to the 
self-interference of the two counter-rotating waves. We model noise as acting on the coin, and given 
by the squeezed generalized amplitude damping (SGAD) channel, which generalizes the generalized 
i-Q ' amplitude damping channel. 
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I. INTRODUCTION 



cr 

Quantum walks (QWs) are the quantum analog of the classical random walks (CRWs) [lMal- It is now well known 
that the quantum features of interference and superposition make QW spread quadratically faster than CRW. Like 
their classical counterpart .QWs are also widely studied in two forms: continuous-time QW (CTQW) [5[ and discrete- 
time QW (DTQW) [3, [j, [a, 0] and are found to be very useful from the perspective of quantum algorithms [81-QjJ . 
Furthermore, they have been used to demonstrate the coherent quantum control over atoms, quantum phase transition 
[l2|, to explain phenomena such as the breakdown of an electric-field driven system [13| and direct experimental 
evidence for wavelike energy transfer within photosynthetic systems [14j . Ila ] . Experimental implementation of QWs 
f—^ ■ has been reported with samples in an nuclear magnetic resonance (NMR) system [161 Il7|. in the continuous tunneling 
f^ | of light fields through waveguide lattices [lj|, in the phase space of trapped ions [1JJ, [20| based on the scheme 
7— I . proposed by [21(, with single neutral atoms in optically trapped atoms Q, and with single photons [2J|. Recently, 
implementation of a QW-based search algorithm in an NMR system has been re por ted [24| . Various other schemes 
have been proposed for their physical realization in different physical systems [25I428J ] . 

Dccoherence in QW and the transition of QW to CRW is quite important from the viewpoint of practical implc- 
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mentation, and has been studied by various authors [29I434J ] . In particular, in Rcfs. J32|. |33|. we investigated some 
qualitatively different ways in which environmental effects suppress quantum superposition in QW on a line and an 
n— cycle. In these works, we have modeled noise as the bit flip channel, the phase flip channel (dccoherence without 
dissipation), and the generalized amplitude damping channel (dccoherence with dissipation), which are of relevance 
to studies in quantum optics and condensed matter systems. The latter type of noise in the QW context appears to 
have been first studied by us [32J. Here we employ a more general noise acting on the coin, the squeezed generalized 
amplitude damping (SGAD) channel, that arises due to coupling with a squeezed, thermal bath [351 ] . 

In this work, we consider the issue of quantifying the degree of quantumness of the walker during the transition 
from a pure QW dynamics to the CRW dynamics. As a general problem applied to an arbitrary quantum system, this 
is still not fully understood. Intuitively, one would expect that quantum entanglement captures the entire content of 
quantum correlations in a system, but this may not be so, in view of other contributions, coming from superposition 
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and noncommutativity of obscrvablcs. In Rcf. [36|, a quantity called 'discord', which is defined as the difference 
between two different quantum generalizations of classical mutual information, is shown to be a viable measure of the 
quantumness of correlations in a bipartite system. One can produce systems in which discord is nonvanishing even 
though entanglement vanishes. However, computing discord in a large-dimensional system can be difficult because 
it requires an optimization over all possible local measurement strategies. An alternative method is via 'dissonance', 
defined as the relative entropy between a separable state a and the closest classical state (a separable state that is a 
mixture of locally distinguishable states), where a is such that it has the smallest relative entropy of entanglement 



with the given (in general, entangled) state p 37j. Therefore, unlike discord, dissonance is a measure of quantumness 
that excludes entanglement. An operational app roach that employs measurement-induced disturbance as a measure 
of quantumness has been proposed by Luo [381 ] . In this work, we use this 'quantumness via measurement-induced 
disturbance' (QMID) between the coin and position degrees of freedom to study the effect of decoherence on QW. 

This article is arranged as follows. In Sec. |TT]we describe the DTQW model on a line and an n— cycle. In Sec. IHII we 
introduce the SGAD channel as the process that decoheres the walk, and discuss some signatures of the classicalization 
of QW. In Sec. HVl we discuss a specific kind of walk symmetry as a means to characterize the loss of quantumness in 
a decoherent walk, and point out its drawbacks. In Sec. |V]we investigate QMID as a means to quantify quantumness 
of the correlations between the coin and the position degrees of freedom, that is applicable to decoherent QWs both 



on a line and on an n— cycle. Finally, we conclude in Sec. IVI1 

II. DISCRETE-TIME QUANTUM WALK ON A LINE AND AN n CYCLE 

The DTQW in one dimension is modeled as a particle consisting of a two-level coin (a qubit) living in the Hilbert 
space H c , spanned by |0) and |1), and a position degree of freedom living in the Hilbert space H p , spanned by \if> x )> 
where i£l, the set of integers. In an n-cycle walk, there are only n allowed positions, and in addition the periodic 
boundary condition \ip x ) = \ip x moc i n ) is imposed. A t-step coined QW is generated by iteratively applying a unitary 
operation W which acts on the Hilbert space T-L c <S> T-L p : 

\%) = W t \* in ), (1) 

where |* in ) = (cos(<$/2)|0) + sin(5/2)e*^|l)) <g> |r/>o) is an arbitrary initial state of the particle and W = U(B <g> 1), 

/ e i? cos(#) e^ sin(6>) \ 
where 1/(2) 3 B = B^^,c, = I -< • ta\ -%i in\ ) (with components denoted Bjf.) is the coin operation. U 

is the controlled-shift operation 

U = |0)(0| ® Y. W* - WVfel + l=L> <!| ® E W* + X )^l- ( 2 ) 

X X 

For an n— cycle U is replaced by 

n—X n—1 

U C = |0)(0| (g) Y, 1^-1 mod n)(^' X \ + |1)(1| ® J2 l^ + ! mod ")^l ( 3 ) 

x =0 X =0 

The probability to find the particle at site x after t steps is given by p(x, t) = (ip x \tr c (\i&!t) (^t\)\ipx) ■ 

III. DECOHERENCE VIA SQUEEZED GENERALIZED AMPLITUDE DAMPING CHANNEL 

A QW implemented on a quantum computer is inevitably affected by errors caused by noise due to the environment. 
Some physically relevant models of noise are as follows: a phase flip channel (which is equivalent to a phase damping 
or purely dephasing channel), a bit flip channel, a generalized amplitude damping channel (T > 0) and an SGAD 
channel, which is used here. Our numerical implementation of these channels evolves the density matrix at each walk 
step employing the Kraus operator representation for them. 

The physical origin and the effect of various types of noise, such as phase flip, bit flip, amplitude and generalized 
amplitude damping, on the DTQW have been studied in detail in Refs. [32|, [33J]. Here we study the behavior of 
decoherent QW subjected to a recently introduced noise process, the SGAD channel [35||, which we briefly discuss 
below for completeness. 



System evolution in the interaction picture, for the case of SGAD, has the following form [30L 1401] : 



j t p s {t) = 10 {N + 1) (a-p'(t)<r+ - ±a + a_p s (t) - ±p s (t)a + <T_ 

-a_a + p s (t) - -p s (t)a_a + 



+ loN U + p s (t)a. 

- -f Q Ma + p s (t)a + - l0 M*<7_p s (t)<7_ 



(4) 



Here 70 = (4w 3 |d| 2 )/(3ftc 3 ) is the spontaneous emission rate, while o+, er_ are the standard raising and lowering 
operators. Also, N = iV th [cosh 2 (r)+sinh 2 (r)]+sinh 2 (r), M = -\ sinh(2r)e l *(2A th + l), where N th = l/(e hu ° /kBT -l) 
is the Planck distribution giving the number of thermal photons at frequency ujq, while r, $ arc the bath squeezing 
parameters 
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The SGAD channel is generated by the master Eq. (U]), with nonvanishing r, <I>,T. The Kraus operators for this 
channel, depicting the system-bath interaction for a time t = A are [35[ 



E = y/¥l 

E 2 = t/V2 



Here, 
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The SGAD channel extends the concept of a generalized amplitude damping (GAD) channel, which arises due to 
the dissipative interaction with a purely thermal bath to that with a squeezed thermal bath. One recovers generalized 
amplitud e dam ping channel, if the squeezing parameter r — > 0. 

Figure 1(a) depicts the gradual classicalization of a QW on a line under the action of the SGAD noise for different 
channel parameters. It is seen that, as the QW turns into a CRW, correspondingly, the probability distribution 
becomes increasingly Gaussian, causing the quadratic functional dependence of variance on time, characteristic of 
quantum behavior, to become linear ( cf. Fig. 2 of Rcf. [301, and the analytical results presented therein). This 
pattern of quantum-to-classical transition has also been observed experimentally, in discrete-time quantum walks 
using single-photon, subjected to decoherence of the pure dephasing type [23|]. This behavior is quite generic, and 
can be shown, for example, to arise under arbitrary Markovian decoherence of a continuous-time quantum walk on a 
grap h (cf. Fig. 3 of Ref. [42| ) ■ The Gaussianization is directly reflected in the fall of standard deviation, depicted in 
Fig. 1(b) because classicalization causes probability to cluster closer to the mean position. As it turns out, neither 
the Gaussianization, nor the concomitant dramatic drop in standard deviation, happens in the case of QW on an 
n-cycle [Fig. 2(a) . Thus the characterization of loss of quantumncss in a QW in terms of shape change of the 



probability distribution according to this pattern is, while convenient for walk on a line, not fruitful in the case of 
cyclic walk. In the following section, we will discuss why a sort of 'gauge freedom' in QWs also is not an appropriate 
geometry-independent measure of classicalization of QW. 
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FIG. 1: The effect of increasing the environmental squeezing parameter of the SGAD channel, on the position probability 
distribution of a particle with the initial state (l/y / 2)(|0) + i|l)) <8> | Vo) and coin toss instruction given by B(0° , 45°, 0°). Other 
channel parameters are fixed at T = 2, 70 = 0.025 and A = 0.1. (a) The probability distribution on a line for different values 
of r, at time t — 100 steps. With increasing noise (parametrized by r), the distribution transforms from the characteristic QW 
twin-peaked distribution to the classical Gaussian, (b) The Gaussianization is evidently accompanied by a reduction in the 
standard deviation, as seen here. The plots are given for different values of the coin toss parameter 6. 



IV. A SYMMETRY OF QUANTUM WALK ON A LINE 

We will call an operation applied at each step of the walk that leaves the position probability distribution p(x,t) 

of the particle unchanged as 'walk symmetry'. The phase gate G(a) — „ ia I acting on the coin is such a walk 

symmetry for walk on line, but not for that on an n— cycle. The physical significance of G is that it helps identify 
a family of QWs that are equivalent from the viewpoint of physical implementation, which can sometimes allow a 
significant practical simplification |32| . 

Suppose the walker begins in the state (a|0) + b\l)) (g> \ipo), with \a\ 2 + \b\ 2 = 1. If the coin operator is denoted by 
the matrix B = {Bjj;}, then after application of this operation followed by the conditional shift, the walker is in the 
state J^ (aBjjO + bBj^lJi, foji—i), where j\ € {0, 1}. Proceeding thus, after t steps the walker on a line and on an 
n-cycle are seen to be, respectively, 



!*£(*)) = E e Uta B JtJt _ 1 ---B njl (B JlQ a + B Jll b)\jt^[2j t -t]), 

31,32,- ■■ ,jt 
I*C(*)> = J2 eUtaB 3 t ,3 t -l--- B 32,3l( B 3l0a + B Jll b)\ ]u ^ [2 j^ t modn] ), 



(5a) 
(5b) 



where J t = ji + J '2 + • • • + 
\ipx) = \ip2J t -t) is given by 



31,32, ■■■ ,3t 

jt and < J t < t. The projection of the state |\&i,(i)), Eq. (|5aj) onto some position state 



(lM*i(*)> = E eiaix+t)/2B h,J t -i ■■■B j2 , jl (B jl0 a + B hl b)\j t ) 

31,32,--- ,3t 

E (e* a(a+t)/2 Bo, J - t _ 1 ■ ■ • B jaJl (B jl0 a + B hl b)\0) 

31,32,--- ,3[t-l] 

+ (e^- 1 ^)/^., ■ • • B htjl (B jl0 a + B hl b)\l] 



where we have used the fact that J t = Jt-i + jt = {x + 1)/2. Therefore the probability to detect the particle at x is 
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FIG. 2: (a) Probability distribution of a particle with the initial state (l/v / 2)(|0) +i|l)) ® |t/>o) and coin toss instruction given 
by S(0°, 45°, 0°) on an n— cycle (n = 51) with increase in value of the r for r = 100. In contrast to the case of QW on a line 
[Figs. 1(a) and 1(b) , note that the classicalized walk does not become Gaussian. All plots have been shifted by the respective 
mean value, and downscaled by a factor equal to the respective standard deviation, in order to bring out the shape change 
in the probability distributions, (b) The breakdown in 'gauge symmetry' is described by nonvanishing Kolmogorov (trace) 
distance between the probability distributions resulting from the nonapplication and application of a phase gate, G(30°) at 
each step of the cyclic walk [33|]. The unitary walk is marked by a persistent, nonvanishing symmetry breakdown. Noise causes 
the gradual suppression of superposition terms, leading to a restoration of the symmetry. 
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independent of a, whence the symmetry. Proceeding similarly in the case of QW on an n-cycle, we find that whereas 
each superposition term accumulates phase in the same way, the spatial spread of the state is restricted by the 
geometry not to exceed n, leading to Eq. (|5b|) . As a result, fixing x fixes Jt mod n = (x + t)/2 mod n, but not J t 
itself, so that the phase terms in the superposition Eq. (|5b|) do not factor out globally in the above manner as with 
QW on a line, whence the breakdown in the symmetry in the cyclic case [Fig. 2(b)] . We use the notation n = 2s + 1, 
where s is the number of vertices on the left or right arm of the cycle, excluding the origin. We define turns r by 
t = ts, where t is walk time or the number of steps. Physically, the breakdown in symmetry can be attributed to 
interference of the clockwise and anticlockwise rotating waves on the cycle. 

In Eq. (|5a[) . note that the proof of symmetry works just the same even if the S's in each operation are distinct 
operators, and, further, even if they arc arbitrary matrices, not necessarily unitary. With this insight, consider any 
CP map represented by the set of Kraus operators {-Ei, • ■ • , E^}- The operators Ej may be unitary, as with a Pauli 
channel like bit flip or phase flip channels, or not, as with the amplitude damping channel or squeezed generalized 
amplitude damping channel. Consider an ?i-fold application of the noise process represented by this CP map on an 
initial pure state p = \^i n )(^m\- This produces 



o'= E E ^---E nP El---El 



(6) 



The state p' can equally well be regarded as a uniform mixture of the d n unnormalized states of the form 
Ei n ■ ■ ■ E^l^in), where a particular sequence Efc = Ei n ■ ■ ■ E^ can be regarded as a 'noise trajectory'. Note that 
the {Eik} are themselves Kraus operators for they can be shown to be positive and to satisfy 
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FIG. 3: (a) Decrease in QMID with increasing squeezing parameter r for a particle walking on a line at time t — 100 steps, 
with the initial state (l/v2)(|0) + i|l)) ® |V>o) and coin toss instruction given by B(0° ,8,0°), where 8 — 30°, 8 = 45°, and 
8 = 60°. A similar decreasing pattern can also be seen in the cyclic case. Note the slight 8 dependence of Q, and the switch in 
this dependence around r = 1.2. (b) Variation of QMID [Q{p)] for QW on an n— cycle (n = 101), with increase in number of 
turns ( r ) of the cycle. With increase in r, Q(p) decreases faster, implying more rapid classicalization. Other environmental 
parameters are 70 = 0.025, temperature (T = 2), and A = 0.1. 



A particular noise trajectory of n steps of a walk, given by (UB) n \^>i n ) could be, for example, 
{UBEi n ){UBEi n _ 1 )- ■ •(JJBEi 2 )(UBEi 1 )\ i f!i n ), where the noise is assumed to act on the coin just before the ap- 
plication of B at each step. This can be rewritten, matrix-multiplying B with the noise operator, as E x |\IJ,- n ) = 
(UB' in )(UB' lni )- ■ ■ (UB'^iUBlJl^in). If the B's in Eq. (J5aj) are replaced by the B n s, we are led by the same 
reasoning to conclude that along each trajectory, p(x,t) is unaffected by the application of G(a). Hence, the average 
overall d n noise trajectories will also leave p{x, t) unaffected. 

We assumed above that the noise acts at each step just before the coin operation, but it is clear that 
it may equally well act just after the coin operation. Furthermore, it may also act after the shift op- 
eration. To sec this we rewrite, using the associative property of quantum operations, a particular noise 
trajectory (E in V ~B){E iri ^U B) ■ ■ ■ (E i2 U 'B)(E h U 'B)|* L ) as E ln {UBE ln _ t ){UB ■ ■ ■ E i2 )(UBE n )UB\^ m ), then as 
E in {UB' in _ 1 )(---UB' ia ){UB' ii )UB\^ L ) and finally as E in (UBE ln _ x )(UB ■ ■ ■ E i2 ){UBE h )\^' L ) and once again the 
above reasoning can be applied to show that symmetry holds even with this model of noise. In earlier works, we 
explicitly verified this for the case when the noise is a phase damping channel [33 ] and a generalized amplitude damp- 
ing channel [33j. Not surprisingly, noise restores symmetry for the cyclic QW, which can be understood simply as the 



result of suppression of quantum superpositions [33 1 



Figure 2(b) depicts the breakdown of the symmetry in noiseless QW on an n-cycle, and its restoration with the 
application of noise. At each time step, symmetry breakdown is quantified as the trace or Kolmogorov distance [43j 
between the two probability distributions obtained by applying or not applying the phase gate on every step. We note 
that while t < 2, the cyclic walk remains symmetric, as the diffusing particle is 'unaware' of its nontrivial topology. 
From Eq. (j5b[) . we see that the beginning of breakdown is when max(2J t — t) > n, which happens when t > n, 
or t > 2 + (1/s). Clearly, even though the symmetry breakdown is potentially a good probe of walk topology, as 
an indicator of quantumness, it is relevant only to cyclic walk, and thus not a geometry-independent indicator of 
classicalization. 



V. QUANTUMNESS USING MEASUREMENT-INDUCED DISTURBANCE 



Finally we consider the classicalization of QW on a line and an n— cycle under the influence of the SGAD channel, 
using QMID as a measure of quantumness of the correlations between the coin and position degrees of freedom. Given 
a bipartite state p living in the Hilbert space Ha <£> Hb, and the reduced density matrices being denoted by pa and 
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FIG. 4: A comparison of the degree of classicalization between QW on a line and that on an n-cycle, via QMID as the measure 
of quantumness. (a) For a fixed setting of the environment-induced decoherence (here T = 0.1, 70 = 0.05, A = 0.1 and 
r — 0) and any given sufficiently large time, the cyclic walk typically becomes classical more rapidly on account of enhanced 
randomization due to the interference of the counter-rotating waves, (b) This is evidenced by the suppression of QMID with 
increase in the number of turns r, keeping time t and all noise parameters fixed. 



Pb, a reasonable measure of total correlations between systems A and B is the mutual information: 

I(p) = S(p A ) + S(p B )-S(p), 



(7) 



where S(-) denotes von Neumann entropy. If pa 
the spectral resolution of the reduced states is 



YIjPa^a an d Pb = ^2jPb^b, then the measurement induced by 



n( P ) = ^n 



n| P iT4 



IT*. 



(8) 
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The state II(p) may be considered classical in the sense that there is a (unique) local measurement strategy, namely 
IT, that leaves H(p) unchanged. This strategy is special in that it produces a classical state in p while keeping the 
reduced states invariant. 

According to Luo |38[, we may consider as a measure of quantum correlations the quantity D[p, H(p)], where D is a 
suitable measure of quantum distance, such as Hilbert-Schmidt distance or Bures measure. If we accept that 7[II(p)] 
is a good measure of classical correlations in p, then one may consider QMID, given by 



Q(p) = I(p) - I[U(p)], 



(9) 



as a reasonable measure of quantum correlation |38[ . 

Our application of QMID to QW to derive the quantumness of coin-position correlations supports this view. We 
numerically determined p in Eq. (|7|). by applying SGAD noise at each walk step, for various channel parameters like 
temperature T, coupling constant 70, A, and environmental squeezing parameter r. The operator IT in Eq. ([5]), and 
thence Q{p) in Eq. ([9]), are obtained by diagonalization of the reduced density matrices pa and ps- 

Figure [3] presents typical plots of QMID in a QW system subjected to a SGAD channel with different variable 
parameters. In Fig. 3(a)| we depict the fall of QMID in noisy QW on a line with increasing environment squeezing r. 
A similar pattern emerges also in the case of cyclic QW. Figure |3(b)| depicts the fall of QMID in a noisy, cyclic QW 
with increasing time, a pattern that is seen also in the case of QW on a line. Both figures demonstrate, as expected, 
that increasing the level of noise (parametrized by r) is detrimental to the degree of quantumness. More generally, 
these findings suggest that QMID provides the desired geometry-independent way to describe the classicalization of 
QW, unlike the previously discussed ways. 

With this result, we may compare the decohering properties of QW on a line with those of cyclic QW (Fig. @|. 
For fixed noise setting, at any given sufficiently large time, as seen from Fig. 4(a)| cyclic walk is more classical 



than walk on the line, which may be understood as follows. Dccohercnce randomizes phase, and in cyclic walk, the 
sclf-intcrfcrence of the two oppositely rotating waves enhances this randomization. This explanation is supported by 
Fig. |4(b)| where increasing the number of turns r keeping time t and all noise parameters fixed is seen to suppress 
the quantumness. 

VI. CONCLUSIONS 

In this work, we investigated the loss of quantumness in a decoherent DTQW on a line and an n-cycle. We showed 
that QMID provides a reasonable geometry- independent description of the degree of classicalization of the walk, unlike 
such familiar criterion as the Gaussianization of the walk probability distribution, which applies only to walk on a line, 
or a specific symmetry principle, which applies only to cyclic walks. The noise is modeled as an SGAD channel acting 
on the coin. One might consider the amount of entanglement between the coin and the position degrees of freedom as 
a measure of quantumness, but, as mentioned in the Sec. HI this might be an underestimate because it ignores the non- 
classicality due to non-commutative observable algebra of the coin and position degrees of freedom. Quantum discord 
is a transparent measure of quantumness, but its computation for large systems like a quantum walker can be hard 
because of the optimization it involves. This suggests that QMID is a suitable, geometry-independent operational tool 
to describe the quantumness of QW. As a first application of this idea, we showed that cyclic QWs tend to classicalize 
faster than QWs on a line on account of more efficient phase randomization. 
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